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Abstract 

(N 

In the framework of Lindblad theory for open quantum systems, we calculate the en- 
■ tropy of a damped quantum harmonic oscillator which is initially in a quasi-free state. The 

r**"- . maximally predictable states are identified as those states producing the minimum entropy 

increase after a long enough time. In general, the states with a squeezing parameter de- 
iy-) ' pending on the environment's diffusion coefficients and friction constant are singled out, but 

if the friction constant is much smaller than the oscillator's frequency, coherent states (or 
^ ! thermalized coherent states) are obtained as the preferred classical states. 
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Although quantum theory has reported during its almost one century of existence spectacular 
successes in explaining phenomena at the atomic scale and below, the relation between quantum 
and classical reality is still a problem not very well understood. At the heart of this problem is the 
superposition principle, which gives counter-intuitive results when applied to macroscopic bodies 
(the infamous Schrodinger cat is probably still the best illustration of this apparent paradox). 

In the last 15-16 years, the theory of decoherence 0, |2| has been advanced as a solution 
to this conundrum: macroscopic bodies are open systems, interacting continuously with their 
environment, and it is exactly this interaction which destroys the annoying superposition terms, 
thus creating superselection sectors in the Hilbert space. 

This strategy has been applied by W. H. Zurek and collaborators @, f|] to study the decoherence 
mechanism in the case of a quantum harmonic oscillator undergoing quantum Brownian motion 
||. In this model, the system is coupled linearly through its coordinate to a heat bath of harmonic 
oscillators which constitute the environment. In their seminal paper Zurek, Habib and Paz 
[ZHP] showed that, if one adopts "predictability sieve" as a criterion for classicality - so that the 
preferred classical states are those leading to a minimum increase in entropy - then coherent states 
are singled out as the maximally predictable pure states. 

Recently, M. R. Gallis || applied ZHP's criterion to a harmonic oscillator weakly coupled 
to an environment; the oscillator is treated as an open quantum system in the framework of 
Lindblad theory [0] . The theory of this system was initially developed by Sandulescu and Scutaru 
H and further continued in ||. One of the advantages of adopting this framework is that the 
Lindblad form for the evolution operator is a rather general phenomenological master equation 
which preserves the positivity of the density matrix and entropy and which encompasses a large 
variety of concrete physical problems, from quantum optics to nuclear physics. 

In this paper we generalize Gallis' result (for Lindblad generators linear in position and momen- 
tum) to the case of an arbitrary large coupling between the system and the environment and show 
that for a damped harmonic oscillator the preferred states have a squeezing parameter depending 
on the environment's friction constant and diffusion coefficients. The calculation is done exactly, 
not perturbatively as in ||, which allow us to consider as a particular case the weak-coupling 
limit in which Gallis worked. We prove that coherent states (or thermalized coherent states, if 
one starts with mixed states) are indeed selected when the friction constant is much smaller than 
the angular velocity of the oscillator, the magnitude of the diffusion coefficients being irrelevant. 

The general form of a completely dissipative mapping L is given by the Lindblad theorem 

Hp) = ~[h, p ] + ^J2({v jP ,v;] + [v^pv;]), (i) 

so the evolution equation for an open quantum system is 

dp(t) 



dt 



L(p(t)). (2) 
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one has to take Vj as a first-degree non-commutative polynomial in coordinate and momentum, 



For a harmonic oscillator, H = f — h ^-q 2 . To get an exact soluble model, as shown in 
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and because p and q span the linear space of the first-degree polynomials in p and q, there exist 
only two linear independent operators V\ and V2, with Vi = ciip + kq, i — 1,2 and at, bi complex 
numbers. With the notations 

h 2 h 2 ft 2 

A« = 2 51 M 2 > d pp = 2 51 N 2 > A*? = Ab> = ~2 Re ^ a * jbil ^ 

j=l j=l j=l 

for the entries of the phase-space diffusion matrix (D pp , D pq and D qq are called diffusion coeffi- 
cients) and 
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A = -/m^a*6 j , (4) 
representing the friction constant, takes the form 

L{p) = -r[H,p]-^[q,pp + Pp\ + ^\p,pq + qp\- 

-|f [p, Ml - |f [g, [q,p]) + 2 -^[p, [q,p]]. (5) 

The evolution equation (JJ) with the dissipative mapping L given by (|5|) has the property that the 
set of quasi-free states is invariant under the action of the evolution mapping p(0) — > p{i). 
Quasi-free states are generalizations of the coherent and squeezed states, and are completely 
characterized by the correlation functions: 

a qq = Tr(pq 2 )-[Tr(pq)] 2 , (6) 
a pp = Tr(pp 2 ) - [Tr(pp)] 2 , (7) 

a pq = Tr(p Pq - qP ) - Tr(pp)Tr(pq). (8) 



Let us denote the correlation matrix by 



o=[ mua qq a (g) 



a — 



Because of the Heisenberg inequality, & qq a pp — a pq = deta > j, cr is a complete positive matrix 
which can be diagonalized [JTTJ 



hA „ T ( K 2 



° = —o 1 ( Q K - 2 ) o, (10) 



where 



A = -VdeUr>l (11) 
n 



is the "area" in phase space measured in units of §; O is an orthogonal symplectic matrix, 

O -- 



cos 9 — sin 9 
sin 9 cos 9 

and K is a positive real number (the squeezing parameter of the state). The degree of purity 
is indicated by A > 1; pure states for example correspond to A = 1. The matrix O describes 
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a rotation in phase space, and K a change of scale (squeezing and dilatation) in the phase-space 
coordinates p and q. Given a quasi-free state, A is uniquely determined, but 9 and K are determined 
only up to the transformations 



8 -> mi + 8. 



[n integer) 



and 



e 



The statistical entropy of a quasi-free state is [T^ 

_ A+ 1 A+ 1 
6 = In 



1 A- 1 
— In . 



(12) 
(13) 

(14) 



The study of the time-evolution of S is equivalent to the study of the evolution of A, i.e. how 
deter (i) changes in time. What we will do next is to calculate deter (t) at t 3> A -1 , and find the 
state <t(0) for which S(t) is minimum (because S increases with A, this will be equivalent with 
deter (t) being minimum). 

The time evolution for the correlations cr pp , er qq and er pq is known (see Eq. (4.55) from |5J) 



with 



(see Eq. (4.51), [§ 
5 is given by 



We have then 



a(t) = K(t)(a(0) - cr{oo))K T {t) + a(oo) 

n(t) = 



(15) 
(16) 



' 5(t) -mu(3(t)\ 
a {t) ) 

\ raw \ ' / 

Now, (4.22) and (4.23) from M imply that the time evolution for a, (3, 7 and 



fa(t) (3(t)\_ p -xt 
\j(t) 6(t)J 



K{t) 



-At 



cos out sin cot 

raw 

mu sin ut cos ut 

cos ut sin ut 
-sin cut cos ut 



(17) 



and using ( |i5|) and the notations 



mu 
cr pp (oo) 
mu 



mucr qq {oo J 



we finally get 

efe£cr(t) = e~ 4A *det[a(0) - (t(oo)] + 
+cfetcr(oo). 



cos 2 ut + muer qq (oo) sin 2 u;£ + 2a pq (oo) sin ut cosc^t 



sin 2 ut + muer qq (00) cos 2 u;£ — 2er pq (co) sin cut cos ut 



cr pp {oo) 



mu 



sin a;t cos ut + er pq (oo) (cos 2 ut — sin 2 wt) 



(19) 
(20) 
(21) 



mu{a qq {0) - a qq (oo))T pp {t) + Tqq (t) - 2{er pq (0) - a M (oo))T M (t) 



my 



+ 
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So, what we have here is an exact expression for the phase-space area "occupied" by the system 
at any time t - and implicitly for the entropy S, via the relations (|TT|) and (|14"D. 

As emphasized by Gallis ||, classicality should be an enduring property and not a transient 
one, so large-time scales should be relevant; in our case, the time scale is set by A -1 . Thus, for 
t ^> A -1 , the term e~ 4Xt det[a(0) — cr(oo)] in the expression of deta{t) can be neglected, and to get 
the maximal predictable states we have to minimize the quantity 



-2\t 



muj(a qq (0) - a qq (^))T pp {t) + M°) a ^°°\ qq ( t ) - 2(a pq (0) - a pq (oo))T pq (t) 



[22) 



Using the decomposition (pTOf) and retaining only the part dependent on the initial state parameters 
in fl2"2"|), we find the expression: 



hA(0) 



{cos 2 0(0) K 2 (0)T pp (t) + K- 2 (0)T qq (t) + sin 2 0(0) K" 2 (0)T pp (f) + K 2 (0)T qq (t) 



-2 sin 0(0) cos 0(0) [*T 2 (0) - K 2 (0)J T pq (t)}, 
which, when minimized with respect to K(0) and 0(0) gives 



tan 20* (0) 



2T 

^^pq 



T -T 



and 



K*(0) 



t x /(t pp - r qq y + 4r p 2 + r pp + r qq 



±J(T -T) 2 + AT 2 +T +T 

-J-y \*pp 1 qq! ~ pq ' PP ' <?9 



1/4 



(23) 



(24) 



(the two values for K*(0) represent the same state, via (|I3| ) and the observation that tan(20) is 



invariant under 9 — > | + 



Now, (Tpp(oo), a qq (oo) and a pq (oo) can be expressed in terms of the diffusion coefficients and 
the friction constant (see equation (3.53) in ||) in the form 



a qq {oo) 



a pp (oo) 



2( m u) 2 X(X 2 + u 2 ) [ ( ^ )2(2A2 + ^ + ^ + WAZ V 
' [ 'mu) 2 u 2 D qq + (2A 2 + uj 2 )D pp - 2mu 2 \D t 



2\(\ 2 + u 2 ) 
1 

2mA(A 2 + u 2 

Inserting these expressions into ( fl~9|) - (pl|) and further into (|24]), we obtain 



'P9 



-X(mu) 2 D qq + XD PP + 2m\ 2 D 



(25) 
(26) 
(27) 



K*(0) 



pq 



A 2 



±\l{muD qq - 2a. 



+ 4DL + Jl 



1/4 



(2? 



This is the general form for the squeezing parameter of the quasi-free states selected via ZHP's 
"predictability sieve". In particular, for A C w, the squeezing parameter becomes K*(0) = 1. 
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This basically means that if the initial states are pure then coherent states will be selected. So, 
the same result as in |4| and || is obtained in the limit of a small friction constant (we do not 
have to restrict the discussion, as in || , to the case in which all the effects of the environment are 
relatively small comparing to the Hamiltonian - the diffusion coefficients can be arbitrary large). 

Another interesting particular case in which coherent states are the preferred states, as can be 
seen from (|28|), is that of an isotropic diffusion in phase space (-^r = muD qq and D pq = 0), a model 
which was previously studied in the framework of Lindblad theory by Ingarden and Kossakowski 
13 . 

These results are also independent of the degree of purity of the initial state. If, instead of 
starting with a pure state, the system is initially in a mixed state {A(0) > 1), then states with 
squeezing parameter given by (|28| ) will be singled out. These states are thermalized squeezed 
states in general, and become thermalized coherent state in the limit A <C to or for = muD„ n 

O 7 mu! qq 

and D pq = 0. 

The author gratefully thanks H. Scutaru for enlightening discussions and invaluable ideas. 
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